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Continuous, coherent sub-terahertz radiation arises when a dc voltage is applied across a stack of 
the many intrinsic Josephson junctions in a Bi2Sr2CaCu20s+5 single crystal. The active junctions 
produce an equal number of I-V characteristic branches. Each branch radiates at a slightly tunable 
frequency obeying the ac Josephson relation. The overall output is broadly tunable and nearly 
independent of heating effects and internal cavity frequencies. Amplification by a surrounding 
external cavity to allow for the development of a useful high-power source is proposed. 

PACS numbers: 07.57.Hm, 74.50.+r, 85.25.Cp 



Continuous, broadly tunable, and coherent sources of 
electromagnetic (EM) radiation are presently unavail- 
able for frequencies within the "terahertz gap", the 0.3- 
10 THz range of crucial importance for many applica- 
tions [l], Q . This gap can be filled by the ac Josephson 
effect intrinsic to a atomic-scale layered superconductor. 
Application of a dc voltage V across a single Joseph- 
son junction between identical superconductors leads to 
an ac Josephson current and EM radiation with the 
same frequency / satisfying the ac Josephson relation, 
/ = fj = (2e/h) V, where e is the electric charge and h 
is Planck's constant Using two-dimensional arrays 

of Josephson junctions between wires of superconduct- 
ing Nb, coherent radiation was observed when the array 
was placed parallel to a Nb ground plane (6| , that ampli- 
fied the radiation, but did not affect its frequency, which 
obeyed the Josephson relation. However, the technical 
problems involved in mass production were formidable. 

The layered, high transition temperature T c supercon- 
ductor, Bi2Sr2CaCu208+5 (Bi-2212), behaves as a stack 
of intrinsic Josephson junctions (IJJs) [7]. In Bi-2212, 
each of the junctions is naturally identical, as they are 
evenly spaced with two junctions per unit cell c-axis 
edge length of 1.533 nm. Recently, continuous, coher- 
ent EM radiation was induced by applying a V across 
the stack of N IJJs present in small mesas milled out 
of single crystalline Bi-2212 (jUli). The mesas typically 
have thicknesses d ~ 1-2 /am and areas that vary from 
~ 4 x 1CT 9 m 2 to ~ 4 x 10 -8 m 2 . The thin mesa shapes 
were mostly rectangular, but some were square or circu- 
lar [l6[. Since Bi-2212 is extremely anisotropic, behaving 
for E || c as an insulator Q, the three-dimensional mesa 
structure also behaves as an internal EM cavity, which 
couples to the non-linear ac Josephson currents generated 
in each junction [8l4l8|. Besides satisfying the ac Joseph- 
son relation for a stack of N IJJs, f = fj = (2e/h) V/N, 
it was also consistently reported that / locked onto an 




0.0 0.5 1.0 0.0 1.0 2.0 0.4 0.5 0.6 

Voltage (V) Output (mV) Frequency (THz) 



FIG. 1: (color online) Rl structure, c-axis resistivity, 
outermost-branch IVC, emission output, and frequency spec- 
tra, (a) SIM image, (b) AFM image along A-B in Fig. 1(a). 
(c) Outer-branch IVC. Inset: R c (T). (d) Plot of I (V ut) cor- 
responding to Fig. 1(c). Open circles: Tt op = Tonset = 84.4 K. 
(e) Spectral intensities from points A-E in Fig. 1(d). 



internal cavity mode frequency where m and p 

are integers appropriate for the geometry. Neglecting 
heating effects, for a very thin (d <C w < £) rect- 
angular cavity of length t and width w, the relevant 
transverse magnetic TM 2 (m,p) modes have frequencies 

fm, P = ( c o/2n) \J (m/w) 2 + (p/£) 2 , where Co is the speed 
of light in vacuum and n = yfe « 4.2 is the index of 
refraction for E \\ c in Bi-2212 [§& 



Most workers have thought the enhancement of the 
output radiation by the excitation of an internal cav- 
ity mode was so strong that the radiation from the ac 
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FIG. 2: (color online) Emission from IVC points at the color- 
coded frequencies was observed at the filled diamonds. No 
radiation was detected at the open diamonds. Arrows indicate 
the numbers N&t of resistive junctions from fits to the ac 
Josephson relation, (a) Rl. (b) R2. Insets: Full IVCs. 



FIG. 3: (color online) Emission frequencies / and intensities 
from Rl, and / versus resistive junction numbers iVfi t for Rl 
and R2. (a) f (V) plots for 11 Rl IVC branches. Dashed 
lines: Fits to / = fj = (2e/h)V/N Rt . Inset: f(V/N ht ) plot 
for all Rl data. Dotted line: / = fj. (b), (c) / (iVfit) ranges 
for Rl and R2, with dotted eye guides, (d) Peak intensity on 
a logarithmic scale versus / for Rl. Calculated f^ iP ranges 
are given at the top. N&t values are color coded. 



Josephson current source alone was too weak to ob- 
serve [§-14, 23-22 1- However, recently the contributions 
to the output power from the ac Josephson current source 
alone and that enhanced by resonance with an internal 
EM cavity source were found to be comparable in mag- 
nitude TH, EH HE Hi |- This created a great deal of con- 
troversy. Here we show clear evidence that the mesas 
can emit radiation at many frequencies, without strong 
interaction with an internal EM cavity mode. More im- 
portant, the resulting radiation is tunable over a broad 
range of frequencies, allowing for the construction of a 
powerful device that could fill the terahertz gap. 

We studied two rectangular mesas, Rl and R2. Rl was 
prepared by milling a groove into the Bi-2212 substrate, 
but R2 was sandwiched between two Au layers. Further 
information on Rl, R2, and three other mesas of various 
shapes are presented in the Supplemental Material jUj. 
Figure 1(a) is a scanning ion microscope (SIM) picture of 
the top of Rl, with w = 99.2-102 /mi, £ = 137-140 /im, 
and groove depth d ~ 1.3 /im, before the electrode at- 
tachments. Rl is ~ 3% longer and wider at the bottom of 
the groove than at the top. Figure 1(b) is an atomic force 
microscope (AFM) scan of the groove profile along the 
A-B line in Fig. 1(a). From this d value, we estimate the 
number of IJJs in the stack to be N max ~ 850. The hys- 
teresis loop indicated by the arrows in Fig. 1(c) is typical 
of the outermost branch of the I-V characteristics (IVCs) 
of an IJJ system. At the largest jump from V = 0.25 to 
0.85 V, all of the iV max junctions switch into the resis- 
tive state simultaneously. In the high-V bias region, Rl 
is inevitably Joule heated at a rate of ~ 16 mW. This 
huge power dissipation may cause the negative differen- 
tial resistance between 1 = 7 and 25 mA, since the c-axis 
quasiparticle resistance R c has a strong T-dependence as 



shown in the inset of Fig. 1(c). On the return region of 
the I-V loop where / ~ 5 mA, we observed several small 
steps, indicating that some of the IJJs make the transi- 
tion from the resistive to the superconducting state. 

Figure 1(d) shows the Rl radiation output intensity 
V ou t generated from / on the same scale as in Fig. 1(c). 
The dashed curve indicates the expected drift due to 
thermal radiation from the sample and its holder. Soon 
after the largest jump in the I-V loop, intense emission 
of EM waves was clearly observed. Strong spatial vari- 
ation of the temperature T t0 p over the top mesa surface 
was previously observed [llf, 13, Ei|, and was thought 



to strongly affect the cavity mode-locking condition. We 
take T top to be the spatial average of the temperature 
on the mesa top. For Rl, T c = 81.1 K with a width 
AT C 6.6 K. Below 

^top — ^onset — 84.4 K, indicated by 
the open circles in Figs. 1(c) and 1(d), Rl makes the tran- 
sition to the superconducting state, as shown in Fig. 1(d). 

Figure 1(e) shows the radiation spectra, offset by 150 
counts for clarity, measured at points A-E in Fig. 1(d), at 
which all of the iV max IJJs are in the resistive state. The 
narrow and intense peaks in the emission spectra of B-E 
have maxima varying from 0.47 THz to 0.53 THz. The 
detector resolution- limited widths of the B, C peaks are 
too narrow to arise from synchronization by an internal 
cavity mode alone, and might involve heating effects [l3| . 
Such outermost branch / tunability of up to 40% was 
found previously by varying both V and the bath T [l3|, 
Til , flit • Although the wide / range strongly violates the 



internal cavity resonance condition, / = the ac 



Josephson relation for a stack of 7V max resistive IJJs , 
f = fj = (2e/h) V/N max , is excellently obeyed. 

We found that EM emission also occurs at many points 
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in the inner region of the multiply-branched IVCs, where 
N = 1, 2, ... , 7V max is fixed but different for each branch. 
The synchronization of the N IJJs on a single I-V branch 
is also not controlled by an internal cavity resonance [26j , 
unlike some predictions [l?} • In Figs. 2(a) and 2(b), the 
radiation frequencies, /, are plotted as color-scaled sym- 
bols on the high-bias regions of the multiply-branched 
I-V structures for Rl at 35.0 K and R2 at 52.5 K, re- 
spectively. The insets show the full IVCs. All of the Rl 
curves in Fig. 2(a) bend backwards with increasing cur- 



rent, indicative of Joule heating [ll|, [13|, LL4|, |28| . However, 
R2 is less susceptible to heating effects, and its IVCs in 
Fig. 2(b) are monotonic. The / spectra were obtained at 
as many (J, V) bias points as possible. At the bias points 
denoted by open diamonds, no emission was detected. 

By repeated measurements of the emission from a par- 
ticular I-V branch with constant A we confirmed that / 
satisfies the ac Josephson relation / = fj = (2e/h) V/N. 
In Figs. 3(a) and 4(a), we replotted the emission data 
from the IVCs of Rl and R2 shown in Figs. 2(a) and 2(b) 
in terms of / (V), representing the data from each branch 
in terms of unique symbols and colors. By fits of the data 
for a particular branch to the ac Josephson relation with 
N = A/fit, the experimental best-fit value Afi t for each 
emitting branch was determined. Three examples each 
for Rl and R2 are respectively indicated by the dashed 
lines in Figs. 3(a) and 4(a) and the arrows in Figs. 2(a) 
and 2(b). In the insets to Figs. 3(a) and 4(a), the entire 
emission data from the inner IVC branches of Rl and R2 
are respectively plotted as fiV/N^t). For both Rl and 
R2, the ac Josephson relation / = fj = (2e/h)V/N&t, 
represented by each dotted line, is well obeyed. 

Since Figs. 3(a) and 4(a) clearly demonstrate that / 
is slightly tunable on each branch (indicated by a fixed 
symbol and color), in Figs. 3(b) and 3(c), we respec- 
tively replotted the data from each Rl and R2 branch 
as / (A/fit). Although for both samples / is tunable as 
indicated by the vertical bars for each fixed branch num- 
ber Afit, Rl and R2 display rather different aspects of 
tunability. The overall tunability of Rl and R2, which 
are respectively tunable from 0.44 THz < f < 0.78 THz 
and from 0.43 THz< / <0.76 THz, is nearly the same. 
However, for Rl, this range is primarily due to the f- 
dependence on I-V branch number, whereas for R2, the 
tunability is greatest on a single branch. 

We measured the spectrum for each emitting point in 
the inner IVC branches of Rl and R2, and determined 
its peak intensity as for points B-E on the outer branch 
of Rl shown in Fig. 1(e). In Figs. 3(d) and 4(b), we 
plotted the peak intensity of the emissions on a loga- 
rithmic scale versus / for Rl and R2, respectively, and 
Afit is approximately coded with color. The ranges due 
to mesa profile variations of the respective internal cav- 
ity mode frequencies f^ p are indicated at the figure 
tops. For Rl, no emission was found for / < 0.40 THz, 
excluding the expected cavity resonance frequencies of 



0.255 and 0.357 THz corresponding to the TM 2 (1,0) and 
T]VP(0,1) modes. Moreover, the spectrum observed in 
Fig. 3(d) is almost completely unrelated to any of the 
internal cavity modes. For example, the strongest in- 
tensity observed for / « 0.575 THz, corresponding to 
V = 0.894 V, I = 9.56 mA, and AT fit = 753, is far 
from the two nearest cavity resonance frequencies. Al- 
though the indicated cavity frequency ranges could be 
shifted to higher frequencies by shorter effective i or d 
values due to hot spots [13| , since T = 35 K in each mea- 
surement, such shifts could still not explain the broad 
range of observed frequencies. For A" max = 851, there 
are ~4.3xl0 130 ways to have N^ t of the A" max junctions 
in the resistive state, each of which could in principle 
lead to a different emission intensity. Hence, multiply 
cycling through the IVCs led to a large variation in peak 
intensity at that (I, V) point. The synchronization of 
the Afit < Amax resistive junctions thus occurs indepen- 
dently of the internal cavity mode excitations [26| . More 
important, the emission spectrum in Fig. 3(d) is contin- 
uous from 0.44 THz < / < 0.78 THz, except for a gap 
between 0.66 THz and 0.73 THz. The gap may be due to 
some experimental difficulty in accessing the appropriate 
branches. This widely continuous spectrum independent 
of internal cavity resonances implies that / is broadly 
tunable. Figure 3(b) suggests that the / values from 
lower branch- number emissions could exceed 1 THz. Fur- 
ther tunability might arise from base T variations [l3| . 
as discussed in the Supplemental Material [25j. 

The peak intensity spectrum of R2 measured at 52.5 K 
shown in Fig. 4(b) is also broadly tunable, with emission 
for 0.43 THz < / < 0.76 THz, and a small gap between 
0.67 and 0.75 THz. For R2, the / range from 0.43 THz 
to ~ 0.55 THz is below the range of the lowest cavity 
resonance frequency /f values. Unlike the peak inten- 
sity spectrum of Rl, the largest intensity at 0.565 THz 
is close to this /f range. In this case, the synchroniza- 
tion of the N resistive junctions could be aided by the 
excitation of this internal cavity mode [26|, [27} . 

We measured the full angular dependence of the radia- 
tion from R2 on the outermost branch at / = 0.639 THz, 
and the results are shown and discussed in the Supple- 
mental Material [25]. Although the observed angular 
dependence is similar to that expected from the low-Q 
tail of the TM 2 (1,0) cavity mode excitation, the fre- 
quency is far from /f , but could be amplified by ex- 
citation of a higher frequenc y m ode, as indicated. We 
thus determined [Til, O, [H, Hif that the angularly in- 
tegrated output power at this frequency is ~ 0.2 /iW, 
somewhat smaller than previously reported (il-[l9|. In 
addition, a much larger emission peak intensity, with an 
estimated overall emission power of 4.8 /iW, was found at 
0.557 THz, which is much closer to /f . Besides the ex- 
citation of the TM 2 (1,0) mode, this strong inner branch 
emission could have been enhanced by R2 being sand- 
wiched between two Au layers [23[. More important, the 
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Hence, it is now relatively straightforward to produce 
broadly tunable, continuous, coherent radiation over the 
range of the terahertz gap. 
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FIG. 4: (color online) Emission frequencies and intensities 
from R2. (a) Plots of / (V) for 8 branches in the R2 IVCs. 
Dashed lines: Fits to / = fj = (2e/h) V/N&t- Inset: Plot of 
f(V/N m ) for all R2 data. Dotted line: / = fj. (b) Peak 
intensity on a logarithmic scale versus / for R2. Calculated 
/to, p ranges are given at the top. N^t values are color coded. 



large range in observed / values well below the j{ range 
is clear evidence that the highly tunable emission itself 
can occur without excitation of an internal cavity mode. 

Although Rl was prepared by forming a groove into 
the Bi-2212 substrate, R2 stood atop a Au substrate. 
The heating effects and EM boundary conditions for Rl 
might be significantly different from those for R2 27]. Al- 
though the synchronization of the A^ t junction emissions 
from R2 is likely to be enhanced by internal cavity res- 
onances [2pj , synchronization of the Nr* emissions from 
Rl could arise from the radiation itself [26| , but is more 
likely to be enhanced by the shunt capacitance in the 
electrical circuit arising from the non-emitting insulating 
junctions in the adjacent Bi-2212 substrate (29| . 

To determine whether the excitation of an internal cav- 
ity mode was an essential feature of the coherent radia- 
tion obtained from mesas of Bi-2212 under the applica- 
tion of V across the stack of IJJs in the mesa, we exam- 
ined the inner branches of the IVCs. We found these in- 
ternal branches to emit radiation over a broad frequency 
range. Hence, we conclude that the primary source of the 
intense, coherent sub-THz radiation is the ac Josephson 
current, and that the internal EM cavity produced by the 
geometrical shape of the emitting mesa is at best of minor 
importance, and for one sample, completely irrelevant. 
Hence, broadband, tunable, continuous coherent radia- 
tion can be obtained in the sub-THz / range from Bi-2212 
mesas. Since the wavelength in vacuum is n ~ 4.2 times 
longer than in the mesa, we propose that a high-power 
device can be constructed by surrounding the mesa with 
an external, tunable high-Q EM cavity 30|, 31]. By ex- 
amining yet lower inner branches, it ought to be possible 
to increase the upper / limit into the 1-10 THz range. 
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